DYNAMICAL SUSCEPTIBILITIES IN STRONG COUPLING APPROACH: 
GENERAL SCHEME AND FALIKOV-KIMBALL MODEL 
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A general scheme to calculate dynamical susceptibilities of strongly correlated electron systems 
within the dynamical mean field theory is developed. Approach is based on an expansion over 
electron hopping around the atomic limit (within the diagrammatic technique for site operators: 
projection and Hubbard ones) in infinite dimensions. As an example, the Falicov-Kimball and sim- 
plified pseudospin-electron models are considered for which an analytical expressions for dynamical 
susceptibilities are obtained. 

Key words: dynamical susceptibilities, dynamical mean field theory, Falicov-Kimball model, 
pseudospin-electron model. 

75.40.Gb, 71.10.Fd, 71.27.+a 



CO 



I 

C 

O 

o 



> 

(N 
OS 

m 
o 

O 

c3 



i 

C 

o 
o 



X 



I. INTRODUCTION 

The development of the Dynamical Mean Field The- 
ory (DMFT), which is exact in the d = oo limit, clarified 
some problems connected with the simultaneous consid- 
eration of the electron hopping and strong local correla- 
tions and stimulated a large progress in the understand- 
ing of the strongly correlated electron systems Q| . It was 
shown by Metzner and Vollhardt QJ^] that in the d — oo 
limit self-energies are single-site quantities (do not de- 
pend on wave vector) which leads to a significant sim- 
plification. But such a local self-energy only probes local 
properties in this limit and cannot detect instabilities as- 
sociated with a specific wave-vector jjj that requires the 
calculation of susceptibilities. 

An analytical schemes developed for the descriptions 
of the strongly correlated electrons can be divided into 
two types: (i) weak-coupling theories, that are based 
on the expansion over local many-electron interactions 
(Dyson approach) and are close to the standard Fermi- 
liquid theory, and (ii) strong coupling theories, that start 
from the expansion over electron hopping around the 
atomic limit (see Ref. Q and references therein). In the 
weak-coupling approach susceptibilities are solutions of 
the Bethe-Salpeter equation and it is known that in the 
d — > oo limit the corresponding irreducible four-vertices 
are also local . 

The aim of this article is to develop a general scheme 
to calculate dynamical susceptibilities within a strong- 
coupling DMFT approach for strongly correlated elec- 
tron systems. A special case of the binary alloy [Falicov- 
Kimball (FK)] type models, that can be solved exactly in 
the limit of infinite dimensions , will be also considered. 
Static susceptibilities for the FK model have been al- 
ready investigated by Brandt and Mielsch || , who found 
an Ising-like phase transition to a chess-board charge- 
density- wave phase at half filling, and Freericks 0], who 
showed that system also displayed the incommensurate 
order and phase separation at other fillings. 



The preliminary short version of this article was 
published in the Proceedings of M2S-HTSC-VI Con- 
ference ||. Here we present the details of the general 
scheme. 



II. TWO APPROACHES IN MANY-ELECTRON 
THEORY 

In general, the Hamiltonian of the electronic system 
with strong local correlations can be written in the fol- 
lowing form 



where 



H — H t + ffioo 



(2.1) 



(2.2) 



describes an intersite electron hopping and H\ oc is a sum 
of the single site Hamiltonians 



(2.3) 



which describe local electron correlations and/or inter- 
action with other local excitations (lattice vibrations, 
(pseudo) spins, etc.). 



A. Weak-coupling approach 



As a rule, the first term in fl2.l| ) is considered as an ini- 
tial (zero-order) Hamiltonian and the second one is con- 
sidered as a perturbation. Such approach is well known 
as the Dyson weak-coupling approach, where the single 
electron Green's functions are determined by the Dyson 
equation 



- — s h —>-©-*- (2.4) 



1 



where 



G a k{u v ) = Gak(^u) + Qak{^u)^ak{^y)G ak {uJy) (2.5) 
1 



IU) V + /J, - tk - Y.akiUv) 



Here 



1 



iu) v + fi — tk 



(2.6) 



and the self-energy Eo-j^cj,,) describes all s catt ering pro- 
cesses originated from the second term in (^lj) and can- 
not be divided into parts by cutting off one zero-order 
Green's function G<,k(u v ) line. 

Also, it is well known that within such a weak-coupling 
approach the electron susceptibilities (charge, spin, etc.) 
can be presented in the form 



(2.7) 



where the type of susceptibility is determined by the end- 
ing parts. The full four vertex is a solution of the Bcthc- 
Salpeter equation 



(2. 



with irreducible four vertex D that, in a same way as 



self-energy, describes scattering processes originated from 
local correlations and cannot be divided into parts by 
cutting off two zero-order Green's function lines. 



B. Strong-coupling approach 

On the other hand, an alternative approach based on 
the expansion over electron hopping around the strong- 
coupling limit can be built PJlCfl. In this case, the single- 
electron Green's functions are determined by the Larkin 
equation PJlO|] 



= + 



(2.9) 



Gaki^u) = £ak(uv) + a ak {w v ) t k Gv k {w v ) (2.10) 

1 



(^) _ ** 

where S^^uv) is an irreducible according to Larkin part 
which cannot be divided into parts by cutting off one 
hopping (wavy) line. 

Now, susceptibilities can be presented in the following 
form 



<WV = 'wv + 'wv — >^vw (2.12) 

is a sum of the chains of hopping lines and full four vertex 
is a solution of the equation 



= 0-EC0 



(2.13) 



and 



are 



In Eqs. ( |2.1l| ) and ( |2.13| ) quantities I 

irreducible verticies which cannot be divided into parts 
by cutting off two wavy (hopping) lines. It should be 
noted that, in contrast to the weak-coupling approach, 
verticies in the strong-coupling approach correspond to 
irreducible many-particle Green's functions. 



C. Connection between approaches 



Prom Eqs. (|J) and (|2~To|) one can get a connection 
between the self-energy and the irreducible vertex func- 
tion according to the Larkin equation: 

^■ak( iuj ") = iuJ v + /i - H ak {iu} v ). (2-14) 

The connection between the four verticies in the weak 
and strong coupling approaches is more complicated. For 
the full four verticies one can get 



(2.15) 



whereas the irreducible four verticies are connected by 
equation 



where D' is a solution of equation 



D' 



= @- \ D \\\ D ' 



(2.16) 



(2-17) 



Here 



-L> -D-^/wO- D 

+ + 

-a- — a — -chwvCt- 

-*- -Da/wD- 



(2.18) 



In the case when the operators A and B are constructed 
only by the electronic operators, then for three and two 
verticies we hnd 




D' 



and 



which can be treated as some kind of sum rules. 



(2.19) 



(2.20) 
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III. DYNAMICAL MEAN FIELD THEORY 

In the case of infinite dimensions d — ► oo one should 
scale the hopping integral according to 



U 



(3.1) 



in order to obtain a finite density-of-states and it was 
shown by Metzner || that in this limit the irreducible 
part becomes local 



or 



Sy CT (T - t') = <%H ct (t - t') 



3<j-fc(Wi/) = E a {u) v ). 



(3.2) 



(3.3) 



Such a site-diagonal function, as it was shown by Brandt 
and Mielsch || , can be calculated by mapping the d — ► oo 
lattice problem (2.1) with intersite hopping (3.1) onto the 
atomic model with single-site hopping 



ff t -> / dr \ dr'Y^CAr- 
Jo Jo 



r')at(r)a CT (r'). (3.4) 



Here ^(r — r') is the auxiliary Kadanoff-Baym field (dy- 
namical mean field) which has to be sclf-consistently 
determined from the condition that the same function 
E a (uj v ) defines the Green's function for the lattice and 
atomic model. The self-consistent set of equations for 
S CT (av) and ( a (uv) (e.g., see Ref. [0 and references 
therein) is the following: 



1 



1 



N ^ E a X {u) v ) - t k Ea 1 (u> v ) - Ca{u> v ) 

= g(»)K,{(„K)}), 



where Gg (u) v , {£ it {uj v )Y) is the Green's function for the 
atomic model (3.4). 

In the same way as it was done by Metzner || for the 
irreducible part E a (u) u ), one can prove that in the d — > oo 
limit all irreducible verticies in the strong-coupling ap- 
proach (Eqs. (2.11) and (2.13)) are also single site quan- 



tities. So, they can be also be calculate d from t he at omic 
model. On the other hand, using Eqs. ( |2.14 )-( 2.2C ) it is 
easy to show that the self-energy and irreducible verticies 
in the weak-coupling approach are also local quantities 
(in Ref. @ that was proved using different approach). 

In order to calculate these single-site irreducible verti- 
cies first of all one has to calculate an irred ucib le many 
particle Green's functions for atomic model (3.4), i.e. 



Q =(Ta\a 2 ala 4 ) a 



(3.6) 



TA a|a 2 



Aq) (Ta[a 2 



and 



n lA — Bl i = (TAqBA - (A 



(3.7) 



(3.8) 



Then irreducible four vertex D for the weak-coupling 



approach can be obtained from equation (2.8), where now 
arrows indicate sin gle-electron Green's functions for the 
atomic model ( |3.5| ) and full four vertex is determined 
from the many-particle Green's function by 



n- 



D 



(3.9) 



for the 



On the other hand, irreducible four vertex 
strong-coupling approach can be obtained from equation 
( |2.13 ), where the full four vertex is determined from the 
many-particle Green's function by 



w\j — rvw 
'wv — k/w. 



or 



(3.10) 



(3.11) 



Now, the thin wavy lines are the dynamical mean field 
C<?(wu) and the thick ones, in contrast to ( 2.12j ), represent 
the sum of chains of wavy lines 



Cx(^v) 



E a (w v ) 



(3.12) 



From the expression for the irreducible many particle 
Green's function (3.7), (3.8) and(3.E) for atomic model 

one can find irreducible verticies 



and 



by equa- 



tions 



and 



^,^. < jwwwv (3.13) 
Xx/w »w\, N K/\/\J K/w 



€H3 



(3.14) 



respectively, that complete the calculation of the irre- 
ducible verticies for the dynamical susceptibilities. 

Finally, for the lattice dynamical susceptibility in the 
weak (2.7) and strong ( p. 11 ) coupling approaches we get 



(3.15) 



Ta[a2/ (TagOzi) + (Tajc^) \Ta 



I3CL2 



and 
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respectively. Here, (. . denote an inverse kernels of 
the corresponding integral equations, 



=GW(w„)GW(^) 1 (3.17) 

— J^^^G cr k{^iy)G a k+q{^i'+m)- (3.18) 



and 



IV. BINARY ALLOY TYPE MODELS 

To test the possibilities of the above approach we con- 
sider a binary alloy type model 



Hi 



PC) 



p" 



p~ 



(4.1) 



where Pf — h±S* for the U = pseudospin-electron 
(PE) model §J, P+ = a, Pr = 1 - c, for a binary alloy, 



and ai 



d}d u P? 



fifi,P- = l- f fi for 



the Falicov-Kimball (FK) model @. The main difference 
between these models is in the way that the averaging 
procedure is performed (a statistical one for the PE and 
FK models and a configurational one for the binary alloy) 
and in the thermodynamical regimes (h = const for the 
PE model, c = const for the binary alloy and n c i = const 
and/or nr = const for the FK model). 

The single-p artic le Green's function for the effective 
atomic model (3.4) is a coherent sum of the Green's 
functions for subspaces S? — ±| and is equal to (see 
Refs. IPi) 



Gi a) K) = 



(P+) 



iu> u +fi- Gx(^) - § 
(P-) 



(4.2) 



+ 



iuj v + (J, — Ccr(^) + 



which allows us to find solutions of Eq. (3.5) for the ir- 
reducible part S a (u> u ) and dynamical mean field („(w„). 

The Fourier transform of the two-electron Green's 
function (Ta\a,2a\a4^j is also a coherent sum of the two- 
electron Green's functions for subspaces 



E 



(4.3) 



[iu v +fi- Co{u v ) T f] [wv> + M - Co'{uw) T §] 



that give s for the irreducible two-electron Green's func- 
tion (g_^) an expression which is a sum of the two sepa- 
rable functions 



->V ™W 



1 \ LOuA-m. CO 



v-\-m wj/i ' -\-rn 



<r> A< T (w 1 /)A (71 (a; I/+m ), (4.4) 



where 



A CT (uv) 



gyfiP^HP- 



(4.5) 



Such a separable form allows us to calculate the inverse 
kernels of the integral equations and to find all quantities 
in an analytical form. 

Within t he w eak-coupling approa ch, starting from the 
expression (4.4) with the use of Eq. (3.9), one can find so- 
lution of the Bethe-Salpeter equation for the irreducible 
four vertex in an analytical form: 



per cr 

<7\ a" 



A £r (w I /) 



(4.6) 



G 2 „(«V')+A //K,//) 



i y 

Aff (ul v )K (J1 {uJ v+ m) 

Gff (uj v )G ai (w v+m ) + Aff (oj v )A ai {uj v+m ) '' 



which must be substituted into Eq. ( J2.8D for the full four 
vertex for lattice and allows to calculate electron suscep- 
tibilities, e.g. for charge susceptibility one can get 



X„ (Wfn ) = OmO 



T — 0(T, q) 



+ AT>™), 



K™(ui m ) = ^y — i . 



where 



(4.7) 
(4.8) 

(4.9) 



Aff(w y )A (T ( 



Xa{u) u W m ) [Aff(w„)Aff(w Iy+m ) - x<j(u) v w m )] 
originates from the last term in ( |4.6| ) and 

A„= (4.10) 

-y- 



Aff(w„)Xff(^0)Xffq(w i ,0) 



/? ^ («„0) + A^ ( w „) [ X ff q ( Wl ,0) - x a (w„0)] ' 



e(T,q) 



(4.11) 



Aff(^) [x„K0)-iM)] 



P X* (w„0) + A£ ( w „) ^(^0) - x a {"M ' 
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XaqfavUm) = ~Jr \^ Gak(Uv)G<rk+q{Uu+m), (4-12) 



V. DISCUSSION 



(4.13) 



But such a diagrammatic weak-coupling approach does 
not allow one to calculate pseudospin and mixed suscep- 
tibilities. 

On the other hand, the strong-coupling approach al- 
lows us to calculate all susceptibilities. To do this, we 
calculate the irreducible many-particle Green's functions 



Tn a (T)al,(T')a,ji(T")) — > 6 m o V, A ff (ay ) (u) v ) 

I irr 

v 

■b« a \K{u v )K{u v + m )+G^\u v )G^\u v+m <\ , (4.14) 



{TS*{T)al{T')a„{T")). t 



(4.15) 



(TS z (r)S z (r')) in . -> S m0 (P+) (P~) , (4.16) 
{Tn a {T)n a ,{T')). m ^ ^A^KJ^A,-^) (4.17) 

l> v 1 

V 

{TS z {T)n a {T')). m ^ (4.18) 
S m oV(P + ) (^>E A -(^)' 



After substituting ( [4.14| )-( p~l"8| ) into Eq. ( |3.16] ) we get 

for the charge susceptibility the same expression (4.7), 
whereas for the pseudospin and mixed susceptibilities 



Xq i^m) = SmO 



— 3m0 



A 2 



T-e(T,q) 

(P + ) (P-) 
T-e(T,q) 



(4.19) 



A S *A„ 
T - 9(T, q) 



where 



A S , = V /(P+) (P- 



(4.20) 



(4.21) 



Expression ( 4.11 ) for Q(T,q) coincide with the one 
obtained by Freericks [Q from the equations for the 
static susceptibilities (uj rn — 0) derived by Brandt and 
Mielsch ||. Irreducible four vertex (4.E) was also derived 
by Freericks and Miller Jl4|] within the Baym-Kadanoff 
formalism and used to find the exact solution for the 
nonresonant Raman scattering for the FK model P^,[jl§ . 



For the bina ry alloy (Faliko v-Ki mball) model (4.1) ex- 
pressions @, (PH) and C20|) define the so-called 
isother mal s usceptibilities jlTfj . Furthermore, the pseu- 
dospin ( 4.19| ) and mixed ( |4.20| ) susceptibilities are only 
static (with factor 5 m o) because the pseudospin opera- 
tor Sf commutes with the Hamiltonian and is an inte- 
gral of mot ion. It should be noted, that the binary alloy 
model (4.1) can be reduced to the Ising type model with 
an effective multisite retarded pseudospin interactions by 
taking trac e of the statistical operator with the Hamil- 
tonian (O) over electron (fermion) variables 0. These 
explains the Ising typ e exp ression obtained for the pseu- 
dospin susceptibility ( 4.19| ), but now the expression for 
the critical temperature Q(T,q) is more compli cate d. 

The expression for the charge susceptibility (4.7) con- 
tains two terms. The first one is static and can be written 
as 



,S*S Z 



(5.1) 



This describes the contribution from the pseudospin sub- 
system to the charge susceptibility. It gives the main con- 
tribution to the static susceptibilities. The second term 
K™ n (uj m ) in (17) gives the pure electron response and de- 
scribes the so-called isolated (Kubo) susceptibility Jl7| . 
In general, terms with the factor S m o give the difference 
between the isothermal and isolated susceptibilities (see 
Appendix in Ref. @). 
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